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Abstract 

The influence of three-particle breakup singularities on s-wave meson produc- 
tion in nucleon-nucleon collisions is studied within the distorted wave Born 
approximation. This study is based on a simple scalar model for the two- 
nucleon interaction and the production mechanism. An algorithm for the 
exact numerical treatment of the inherent three-body cuts, together with its 
straightforward implementation is presented. It is also shown that two often- 
used approximations to avoid the calculation of the three-body breakup are 
not justified. The possible impact on pion production observables is discussed. 
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I. INTRODUCTION 



Interest in studies of pion production in nucleon-nucleon collisions at energies near the 
pion production threshold has been revitalized by the appearance of excellent high quality 
data in . In addition, experimental data for double polarization observables for the reactions 
pp ppn^ 0) pn'K'^ 1^, and pp — > dTT~^ ^ are available. A comparison of the 

experimental information with the predictions of a modern meson exchange model shows 
that the data for the production of charged pions are well reproduced, whereas there are 
large discrepancies between the predictions and the data when considering neutral pion 
production. Thus, although the process of pion production has been investigated over several 
decades, the phenomenology of this fundamental process is still not fully resolved. 

A priori one could assume that the most rigorous approach to pion production in nucleon- 
nucleon collisions is a complete description of the nucleon-nucleon-pion system, e.g. via 
coupled channel equations, as developed recently in Ref. [Q. However, not only is such a 
procedure difficult to implement, arguments based on effective field theories suggest that, 
at least close to the production threshold, most of the diagrams generated within a coupled 
channel approach are suppressed Those arguments indicate that the only piece that 

needs to be treated non-perturbatively is the nucleon-nucleon (NN) interaction, whereas 
the transition amplitude NN — > NNtt can be treated perturbatively. In this spirit, most of 
the recent calculations for pion production in NN collisions are based on the distorted wave 
Born approximation (DWBA) Q 

Above the pion production threshold, where a virtual pion exchange can produce a real 
pion, one encounters a final state containing three real particles and therefore has to consider 
a full three-body breakup amplitude. This amplitude contains branch cut singularities. Since 
the evaluation of those so-called three-body singularities is technically rather involved, they 
have been avoided so far in all calculations that apply DWBA to meson production in A^A^ 



collisions [|TT|. Instead, the energy-dependence of the pion propagator was manipulated in 
such a way that no singularity occurred. These approximations, however, were never tested 
quantitatively. Moreover, since the imaginary part arising from the three-body singularities 
scales with the three-body phase space, a strongly energy-dependent phase motion of all 
those amplitudes with a significant contribution from pion exchange should be expected. 
One might speculate that the insufficient treatment of the three-body character of the pion 
production may be a source of the insufficient description of the polarization data for neutral 
pion production mentioned before. 

Recently a simple, scalar 'toy' model for pion production was proposed in order to in- 
vestigate the validity of several approximations used with respect to the energy dependence 
in the pion propagator as well as in the ttN amplitude ||T^. Although the model used in 



Ref. [|12| is simple with respect to the interactions, the pion dynamics is treated exactly. 
A significant finding of that work was that a proper treatment of the dynamics is very im- 
portant, even exactly at the pion production threshold. Therefore, this model is a natural 



^Note: Even before using chiral perturbation theory, DWBA was used for calculations of the 
process A^A^ — > NNtt In these calculations the only justification for this procedure was based 
on the slow growth of the A^A^ inelasticities. 
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choice for studying the impact of an exact treatment of the three-body dynamics in the 
pion production process. In this work we specifically want to address the role of three-body 
singularities in reactions of the type NN NNtt, and to study their possible effect on 
observables. In addition, we present an algorithm for their evaluation that can be imple- 
mented quite easily. For this specific study an exact treatment of the A^A^ interaction is not 
necessary. The inclusion of a full A^A^ interaction will add to the computational complexity 
by requiring careful interpolations of the half off-shell A^A^ transition amplitudes, but will 
have no impact on the technical peculiarities induced by the three-particle singularities of 
the pion propagator. Such a model is thus the ideal starting point for our investigation. 

This article is structured as follows: In Section II we briefly summarize the most impor- 
tant features of the model from Ref. |jl2|; in Section III the algorithm for the exact numerical 
evaluation of the three-body breakup singularity is given; in Section IV we present our nu- 
merical results; and we conclude in Section V. 



II. THE 'TOY' MODEL FOR PION PRODUCTION IN NN COLLISIONS 

Here the most important features of the scalar model for the production of pions in A^A^ 
collisions as introduced in Ref. [|T^ are summarized and discussed. The main features of 
this model are: 

i) The pion is represented by a scalar-isoscalar Klein-Gordon field that couples via a 
Yukawa coupling to the nucleons. Only the emission of this particle is considered. 

ii) The nucleons are treated as distinguishable. As a consequence, pion emission needs 
to be considered only from one nucleon. The symmetric term wherein the pion is emitted 
from the other nucleon is omitted. The simplicity of the model is retained by allowing the 
pion to couple to only one of the nucleons. Therefore pion exchange between two nucleons 
is not contained in the model. 

iii) Only pion rescattering by one nucleon is considered. This pion rescattering is de- 
scribed by a vrA^ seagull vertex which is inspired by the chiral vrA^ interaction Lagrangian. 

iv) The nuclear interaction is modeled by the exchange of a heavy scalar-isoscalar Klein- 
Gordon field, named a, which couples to the nucleons via Yukawa coupling. Since the 
strength of this coupling does not influence the three-body dynamics of the pion production 
process, only small couplings are considered for the one-cx exchange. 

It is typical to treat pion production near threshold using a non-relativistic expansion 
in the nucleon momenta. In the following only leading terms in this expansion will be 
considered; in particular, contributions from anti-nucleons are ignored. 

In summary, we consider a scalar model defined by the following Lagrangian: 

+^4iV27r + g.Y. N]^^^ + y^NlN^{doT^f . (2.1) 

Here M^v represents the physical nucleon mass of 939 MeV and m,r is 139 MeV. The mass 
of the a meson is chosen to be 550 MeV. All diagrams are evaluated at order j-g^^. In 
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the following we do not display these factors, nor other constants that are common to all 
amplitudes. All the loops are finite, so that we do not need to introduce any regulators. 

It is important to point out some of the unrealistic features of the model introduced 
above, even so they will not influence our study of the three-body dynamics of the pion 
production process. We are concerned with near-threshold kinematics, so that a scalar 
particle is produced in an s-wave, as is the final A^A^ pair. Angular momentum conservation 
requires that the initial A^A^ pair is also in an s-wave. In a realistic case, however, where 
the pions are pseudoscalar mesons, the production of s-wave pions calls for a p-wave in the 
initial state. Furthermore, the scalar model does not include a strong short-range repulsion. 
Thus the nucleons have stronger overlap within our model, as compared to a more realistic 
treatment. However, since the main focus of this work is the three-body dynamics of the 
production operator, a study within this model will still give valuable information. As an 
aside, the consideration of more realistic interactions would not introduce new conceptual 
aspects to the exact treatment of three-body singularities, but only increase the numerical 
effort. 

In a DWBA calculation of threshold pion production, any tree level diagram is modified 
substantially by the contributions from the initial and final state interactions. We will 
therefore concentrate in our model calculations on those DWBA terms, in which only initial 
or final state A^A^ interactions are present. We will thus ignore the rescattering diagram 
with DWBA contributions in both initial and final A^A^ interactions, since this is a two-loop 
integration term. Again, the initial and final state interactions, which occur before or after 
the pion rescattering process, are represented by a single a exchange. All diagrams are 
evaluated in time-ordered perturbation theory (TOPT). The diagrams involving an initial 
state interaction are depicted in Fig. In the present work we will only study ladder 
diagrams, which are given by the contributions labeled II and 12 in that figure. There are 
two additional types of diagrams, namely the stretched boxes (Fig. 0-13, 14), and graphs in 
which a cr is exchanged in between the emission and rescattering of the virtual pion. We 
ignore those here. Although both groups have three-body singularities, they introduce no 
additional complications and are thus not relevant for the present study. 



III. CALCULATION OF THE THREE-BODY BREAKUP 

As mentioned at the end of the former section we will study only ladder diagrams. The 
corresponding diagrams are shown in Fig. 2 together with the relevant momenta. In time- 
ordered perturbation theory the corresponding amplitudes are proportional to 

A^ = l dV Vf Gf^^iC:';,^ - Glt^^) , (3.1) 
= j dY {G%^ - Gitr^MNV^ =■■ A^'^ + A^'" . (3.2) 

The superscript r (6) denotes the diagrams for pion rescattering (backscattering) , namely 
those with a pion moving forward (backward) in time, as shown in Fig. 1 II (12). The 
individual pieces can be directly deduced from the diagram; e.g. the two propagators in 
Eq. ( p.2| ) are given by 
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^nNN — \ ^/ T^, ~~) ) (3-3) 

\Etot --£/+- E" -Ug + ieJ 

and 

^''^^ = {E,^t-E'_-E^-u^,-u,,+^e) ' ^^'^^ 

The relative minus sign between the two propagators in Eqs. ( p.l| ) and ( |3.2D stems from 
the particular form of the irN —>■ nN transition operator, namely the pion energies that 
appear explicitly. The evaluation of the expressions for as well as A^'^ is straightforward 
since they only contain two nucleon singularities, as long as total energies Etot below the 
two-pion threshold are considered (c.f. Eq. ( |3.4| )). Thus in the following we will concentrate 
on the evaluation of the amplitude A^'^. When discussing the results, however, the complete 
amplitudes A^ and as A^ will be considered. 

For completeness we will also give explicitly the expressions for the other parts of A^'^, 
namely the two-nucleon propagator 



and the a-exchange potential 



uj^{Etot-E-E''-u^)- ^^-^^ 

In a model based on a realistic NN interaction, the potential will have to be replaced 
by a A^A^ t-matrix. 

The integration variable p" in Eq. ( p.2| ) stands for the relative momentum of the nucleons 
in the intermediate state. The quantity E'^ = jM^iP'^^y energy of the right (+) and 



the left (— ) nucleon in the final state, Uq' = Jm^ + g'^ is the energy of the produced pion 



and the total energy is given by Etot = 2-E = E^^, + E'_ + ujqi. In addition, ujq = \Jm1 + 

I — * -'112 

and uOfj = \J rn^ + kP' denote the on-shell energies of the vr and a meson, and E" = is 

the energy of a nucleon in an intermediate state. The momenta q and k are defined through 
q = p' + \ — p" and k = p — p" . The momentum of the initial nucleon is labeled p (c.f. 
Fig. I). 



The amplitude given by Eq. ( |3.2| ) exhibits a distinct singularity structure; it contains 
both, a two- and a three-body cut. Those cuts occur whenever it is energetically allowed for 
two or three particles in the intermediate state to go on-shell simultaneously. It should be 
noted that the two and three-body singularities are kinematically well separated, as shown 
in detail in Appendix A; whenever the intermediate two-nucleon state is on-shell, the nNN 
propagator is strictly negative. The numerical treatment of two-body singularities, e.g. by 
subtraction methods, is a well-established procedure [jl3[. Calculations of the three-body 
breakup have long history in neutron-deuteron scattering (e.g. |14- [l^] ). In the context of 
NN potentials that incorporate pion production explicitly, the breakup has been treated 
either by complex contour deformation Ijl^, spline methods ||18[, or subtraction methods 
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19| , |20[] . For the explicit calculation of pion production observables the last seem preferable, 



since all calculations are carried out along the real axis. For our purpose, namely the 
framework of DWBA, the method introduced by M. Schwamb pO|j2T[] is most suited. Thus, 
this algorithm is presented here together with modifications necessary for the evaluation of 
the loop given in Fig. 2. 

In order to proceed, the ttNN propagator introduced in Eq. (|3.3| ) is rewritten as 

^I,r _ ( Etot- E'^- E" + UJq \ I 

— » -*/ 

Here P' = p' + ^ cind 



'ml + P'' +p"^ -2P'p"x" , (3.8) 
where x" = {P' ■ p")/{P'p"). The quantity Xq is defined as 

_ ml + P"+p"'-iE,,,-E'^-Ey 
xo = xo{P ,p ) = ^j^, . (3.9) 

The integration over the angle variable x " contains a singularity for all p " that leads to 
|xo| < 1. Thus, 

. ±1 « „5 + (p. „T - (3.10) 

represent two conditions for the boundaries of the area containing three-particle singularities. 
These are two polynomial equations of degree four for p" = p"{P'). Only three of the eight 
solutions are physically relevant, as can be easily checked by taking the threshold limit for 
the solutions. They enclose the area of three-particle singularities, as illustrated in Fig. |^, 
for a total energy of Et^t = 210 MeV. In what follows we will denote the boundaries of this 
area by pa and pb, respectively, as indicated in the figure. 
Introducing the auxiliary function 

F{x") := {E,,t ~E'^- E" + u,) / d^" (3.11) 

allows to rewrite the amplitude A^'"^ from Eq. ( |3.2|) as 

A^'^ = A{^^^ + y4^Q + , (3-12) 

where 

IP' ^ p^-p"^ + ieJ-i x"-Xo + ie ^ ^ 

The regions of integration are given by 

R3b=[Pa,Pb\; Rno=[0,Pa] U [pb,P'i\; R2b = [Pi, Oo] . (3.14) 
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Here pi denotes an in-principle arbitrary momentum placed in the interval between the 
region of three-body singularities and the two-nucleon singularity (c.f. Fig. |^). For the 
numerical evaluation we chose Pi such that it cuts this interval in half. 

The term is thus the piece of the amplitude that contains the three-body cut. The 
second term in Eq. ( |3.12| ) stands for the remaining integration from pi to infinity and 
contains only the two-nucleon singularity. Since this calculation involves only a standard, 
single subtraction it will not be discussed in detail. The third term, A^^, contains the 
integration outside the region of three-body singularities for momenta smaller than pi and 
will be discussed in detail at the end of this section. 

We begin our discussion with the term Ag^. For each value of Xq we subtract and add 
the value of at the singularity and get 



^3b 



M, 



TV 



Pb 



dp"- 



P 



Pa P^ -P 

1 ^ „F{x")-F{xo) 



dx 



X 



Xo 



(^In 


1 - Xo 


— irc^ 




1 + Xo 





(3.15) 



The complication specific to three-body singularities is the appearance of removable loga- 
rithmic singularities: the arguments of the logarithm vanish at the integration limits pa and 
Pb, where |xo| = 1. 

The last piece of Eq. ( p.l5| ) is the imaginary part mentioned in the introduction. It leads 
to an additional imaginary contribution, which is necessarily energy-dependent. 

Let us concentrate now on the logarithm in the last term of Eq. ( |3.15| ) . Here we employ 



a transformation first introduced by Schwamb |^T| that allows to carry out the integrals 
numerically in a quite straightforward fashion. In order to apply the transformation the 
range of integration over p" is split into two pieces, so that there is only one singularity in 
each of the two intervals. The intermediate momentum is chosen to be pab '■= |(Pa +Pb), in 
accordance with Ref. ||2l|. Using the abbreviations a{p") and P{p"), we are now faced with 
integrals of the type 



Vb 



dp"a{p") HP{p")] 



Pab 



where /9(pb) = 0. The substitution 



2 " \ n 

y ■.= pb-p >0 



allows to transform the integral of Eq. ( |3.16| ) into 



I 



Pb-Pa 
2 



dyya{pb - y'^) ln[/3(pb - y"^)] 



(3.16) 



(3.17) 



(3.18) 



Due to the minus sign in Eq. ( |3.17| ) the lower integration limit is now the 'critical' point 
of the integrand. The additional factor y in Eq. ( |3.18| ) is responsible for the 'disappearance' 
of the logarithmic singularity; the integral is now regular, since \im.y^Qy\n{ii) = 0. The 
integration interval [pa,Pab\ is treated in a similar fashion, as described in Appendix B. 

In principle, the application of a subtraction in the intervals [0,Pa] and is not 

necessary, since singularities occur at the boundaries only. However, numerical tests showed 
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that applying a subtraction in these intervals also allows a smaller number of grid points in 
obtaining a stable and convergent result. In summary, with the above algebraic transforma- 
tions, the second term in Eq. ( |3.12| ) is expressed as 



AL 



N 



2P' 



dp"- 



P' 



dx 



,F{x") — F(signa;o) 



X' 



Xo 



+ F(signxo) In 



1 - Xo 



l+Xo 



(3.19) 



where the integration interval Rno was defined in Eq. ( p.l4| ). 

Note that, through the splitting into intervals of the p" integration the first term in 
Eqs. ( p.l5| ) and ( p.l9| ) also appears as integrated in intervals. In principle there are no 
numerical difficulties with this integral, however one needs to make sure that there is a 
sufficient number of mesh points within the region of the three-body singularities in order 
to obtain accurate results. That is why we propose to use the splitting of the intervals here 
as well. A complete presentation of the three-body cut algorithm is given in Appendix B. 



IV. EFFECT OF THE THREE-BODY BREAKUP ON THE PRODUCTION 

CROSS SECTION 

Since the ingredients of our model are of scalar nature, the only observables to be con- 
sidered here are cross sections. The total cross section atot for pion production is calculated 
from the coherent sum of the production amplitudes and A^ defined in Eqs. ( p.l|) and 
( p. 21) , respectively. Here we will use the notation 

a,,, = Jdp\A^ + A^\' ; af^, = J dp\A^\' ; al, = Jdp\A^\' , (4.1) 

where dp denotes the phase space factor. As a measure of the total energy we introduce 
the excess energy Q = Etot — ^tt- In the present work we concentrate on excess energies 
between the one-pion and two-pion threshold. Thus, if the excess energy is given in units of 
m^r, we restrict ourselves to values of Q/m.^ between and 1. In the upper panel of Fig. ^ 
we display the result of the full calculation, atot (solid line), with cross sections (dashed 
line) and cr^^ (dotted line) obtained from the initial and final state interactions individually. 

In order to emphasize better the contributions of the FSI and ISI separately, we divide 
al^^ and cr^^ by atot and display those ratios in the lower panel of Fig. |. The figure shows 
that al^^ gives the largest contribution over the entire energy interval under consideration, 
while cr^^j is always close to atot- 

Naturally these features are specific to the model. More interesting is to compare the 
exact calculation with approximations for the pion propagator to avoid the appearance of 
the three-body singularities, since a main goal of this work is to investigate the quality of 
those approximations. Here we consider two of these approximations. In both, the nNN 
propagator Eq. (|3.3| ) is substituted by a propagator containing no three-body singularity. 
In case of the static approximation (used, e.g., in Ref. [^) the substitution 



Etot - -E+ - E" -ujg + ie Uq 
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is made; i.e., the vr-exchange is made instantaneous. In case of the ^threshold approximation^ 
(used in Refs. the vrA^A^ propagator is given by 

Etot -K- E" ~ujg + te ^ m^-E'+- E" - ujg ' ^^'^^ 

By definition, exactly at threshold this 'threshold approximation' agrees with the exact prop- 
agator, since the total energy Etot is then equal to m^r. In the literature other approximations 
can be found. For example, in Ref. [0 a propagator similar to that of Eq. ( |4.^ ) is used, but 
with E'j^ = 171^^/2 and E" = 0. Recognizing that our list is not complete, we nevertheless 
consider the 'static approximation' and the 'threshold approximation' as representative and 
concentrate on these. 

In Fig. ^ the ratios of the calculations based on the two approximations and the exactly 
calculated cross sections are displayed. In panel (c) of Fig. ^the results of the full calculations 
(i.e. containing FSI and ISI) are shown. The dashed line gives the result for the 'static 
approximation', while the dotted line represents the 'threshold approximation'. The figure 
indicates that, although the 'threshold approximation' is by definition exact at threshold, its 
energy behavior does not correspond to the exact calculation at higher energies. Especially 
when approaching the two-pion threshold, this approximation underestimates the full result 
considerably. The 'static approximation' overestimates the full result in the entire energy 
regime between the one-pion and two-pion thresholds. This is especially pronounced near 
the one-pion threshold. In order to investigate the effects of the approximations in more 
detail we display their effects calculated with FSI and ISI only in panels (a) and (b) of Fig. |^. 
Comparing the panels (a) and (b) with the complete calculation (c) shows that the last is 
the result of interference effects: due to their singularity structure and have different 
phases with different energy-dependence that are taken into account in different ways by 
the two approximations. 

In order to study this point in more detail let us now consider the phase motion intro- 
duced by the exact inclusion of the three-body singularities. Obviously, in order to obtain a 
quantitative understanding, we need to compare the full result to one in which the breakup 
does not occur. For this comparison we choose the 'static approximation', defined in Eq. 
( [4.2| ). In Fig. ^ we compare the phase (plA] of the amplitudes, defined through 

for a fixed total excess energy Q = ^rrij^ as a function of the kinetic energy of the relative 
motion of the outgoing two nucleon system. The solid line shows the phase of the full 
result, while the dashed line corresponds to the 'static approximation', panel (a) contains 
the phase panel (b) the phase and the phase of the full calculation, 0[y4-^ -|- A-'^] 

, is shown in panel (c). As could have been expected, the phase of A^ is infiuenced most 
by the approximation, and there is a significant effect on the phase of the total amplitude. 
We would like to point out that for each value of Q that is smaller than m7r, the effect is 
qualitatively similar to the one shown in Fig. ^. Naturally, this effect is irrelevant as long as 
only one partial wave contributes. However, as soon as differential observables are analyzed, 
different partial waves start to interfere. It will be important to see how such a change in the 
phases of some amplitudes influences the description of the polarization data of the reaction 
NN NNtt. 
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V. SUMMARY AND CONCLUSIONS 



The effect of tliree-body dynamics in tlie pion production process witliin tlie framework 
of the DWBA approximation has been investigated using a scalar 'toy' model. This model is 
simple with respect to the interactions employed: the NN interaction is represented by the 
exchange of a scalar cr-meson; the pion is also treated as a scalar, and finally the nucleons 
(represented by scalar fields) are assumed to be distinguishable. However, the imderlying 
dynamics is treated exactly, especially the three-body breakup NN — > NNt:. We presented 
an algorithm for evaluating the three-body breakup singularities exactly. This algorithm 
can be applied to realistic NN interactions in a straightforward way. 

Within our model we compared the results for the total production cross section with two 
often-used approximations: the 'static approximation', which makes the pion exchange in- 
stantaneous, and the 'threshold approximation', which fixes the total energy at the threshold 
value. It turns out that both approximations are different in character, but equally unsuited 
to describe the total pion production cross section between the one- and two-pion thresholds, 
at least within this model study. In case of the static approximation the cross section is 
overestimated over the entire energy regime, but especially close to threshold. Although the 
threshold approximation is exact at the threshold, it overestimates the cross section slightly 
close to threshold and underestimates it considerably near the two-pion threshold. Both 
approximations produce the wrong energy-dependence in the energy region under consider- 
ation. It remains to be seen how much of these differences survive in a realistic calculation, 
where, e.g., the final state interaction plays an a lot more prominent role and small momenta 
are suppressed due to chiral symmetry. 

Finally, we compared the phases of the amplitudes of our exact calculation with the 
phases given by the 'static approximation' and found sizable differences. Although not 
presented in this work, a quite similar statement can be made with respect to the 'threshold 
approximation'. Naturally, the phases play an important role when different partial waves 
interfere, as occurs for differential observables. It can then be expected that in the case of 
realistic interactions the description of some observables will depend crucially on the way 
that the three-body dynamics of the intermediate states is treated. 
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APPENDIX A: LOCATION OF THE TWO AND THREE-BODY 

SINGULARITIES 



We demonstrate here that the two- and three-body singularities from Eq. (|3.2| ) are well 
separated. The two particle singularity occurs at 

p" = p= ^M^Etot, (Al) 

i.e., in case of on-shell scattering. An estimate for the denominator of the three-body 
propagator clarifies the position of its singularities, namely 

V p>p: Etot- E'^- E"{p" =p) - ujq{p" =p) 

<^-E'^-uj,{p"=p) < 0, (A2) 

since E"{p" = p) > E = and Etot < 2m^, ojq > 171^^, and E'^ = > 0. Consequently 
the TrA^A^-propagator is regular for all p > p. In particular, the position of the two-nucleon 
pole occurs above each three-particle singularity. This leads to the natural decomposition 
of the radial integration interval /o°° = Jq' + (c.f. Eq. ( p.l2|) ), where [0,pi] contains all 
three-body singularities and [pi, 00] contains the two-body singularity. 



APPENDIX B: ANALYTIC AND NUMERIC TREATMENT OF THE 
LOGARITHMIC SINGULARITIES 



Here we give a brief but complete presentation of the three-body cut algorithm. With 



xo if |xo| < 1 
signa;o if \xo\ > 1 



(Bl) 



(i.e., ^0 = xo in Rsb = [pa,Pb] and ^0 = signa;o in Rno = [0,Pa] U [pb,Pi],) the contribution 
+ (c.f. Eqs. ( p.l2|) , ( p.l5|) , and ( p.l9|) ) to the production amphtude reads 

(B2) 



where 



lm{Ai 



7tM 



N 



3bJ 



2P' 



P" „p"F{xo) 
dp — Til ' 

Pa P'^ -p" 



A 



IM 



N 



2P' 



Pa 



+ 



Pab 



+ 



Pa 



Pb 



+ 



Pab 



Pb 



dp 



0) 



In 



Xo 



1 + Xo 



and 



(B3) 
(B4) 

(B5) 

p^ — p"" ■ ' 

^^,in (;ienotes the real contribution to + containing the logarithm originating in the 
principal value of the analytically solvable integral, and A^''^ is the other real contribution 
originating in the application of the subtraction method. The expression A^'^^ =: Ipf Z]j=i 



A 



/, A 



M, 



N 



2P' 



Pa 



+ 



Pb 



+ 



Pa 



Pb 



dp"- 



P' 



,//2 



-1 x" — Xo 



from Eq. ( p4D is to be understood as if the corresponding Schwamb's substitution would 
have been already performed in each of the four integrals (c.f. Table |). 
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TABLES 



TABLE I. Schwamb's substitutions in the four parts of the interval [0,pj]. 



interval transformation contribution 



[0,Pa] 


y' 


■=Pa 


-P" 


> 


-2j^dyyaipa - y^) ln[/3(p, - y^)] 


[Pa,Pab] 


y' 


:=p" 


-Pa 


> 


2 j^^-^'-P^ dy ya{pa + y^) Hf3{pa + y^)\ 

-'^fvKiPp-Jyy^ipb - y") Mi^^^ - y')\ 


\Pab,Pb] 


y' 


■=Pb ■ 


-p" 


> 


\Pb,Pi] 


y' 


■.= p" 


-Pb 


> 


2 /o^P^-^" dy ya{p, + y^) ln[/3(p;, + y^)] 
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FIGURES 



FIG. 1. The six TOFT diagrams of the initial state interaction. A pion field is represented by 
a dashed line, the nucleon by a single solid line. The two possible time orderings of the cr-exchange 
(solid double line) lead to the identical contributions grouped in II and 12, for each of two ladder 
diagrams. 13 and 14 are stretched boxes. 

FIG. 2. Feynman diagram and choice of coordinates for a) final- and b) initial-state interaction. 
All momenta indicated are three-momenta. 

FIG. 3. The region of three-body singularities for £^^(=210 MeV and the choice of the integra- 
tion intervals for the three-body cut algorithm. Inside the region enclosed by the solid line we have 
|xo| < 1 and on the boundary |xo| = 1- Additionally, the two-particle singularity at p ~ 444 MeV 
is marked. The symbol -|- (— ) means that the corresponding boundary curve is solution of xq = 
(xo = -1) (c.f. Eq. (|3l0|)). 

FIG. 4. The total production cross section as function of the excess energy, which is given 
in units of m,r- In the upper panel (a) the full calculation is represented by the solid line, the 
calculation with the FSI (ISI) only is shown as dotted (dashed) line. The lower panel (b) shows the 
ratios of the calculations containing the FSI (dashed) and ISI (dotted) only to the full calculation 
(solid line in panel (a)). 



FIG. 5. The ratios of the total cross section of the static approximation (dashed curve) and 
the 'threshold approximation' (dotted curve) to the total cross section of the corresponding exact 
calculation. The calculations of panel (a) contain only the FSI, the calculations of panel (b) only 
the ISI. In panel (c) the ratios of the full calculations are shown. 

FIG. 6. The phase of the production amplitude in degrees as a function of the kinetic 
energy of the relative motion of the A^A^-system (e in units of mjr) at the excess energy 
Q = ^max = \ fn-K — 70MeV. The solid lines represent the exact calculation, the dashed lines 
the 'static approximation'. In panel (a) calculations with the FSI only are shown, in panel (b) 
those with the ISI only. The full calculations, i.e. the superposition of (a) and (b), are shown in 
panel (c). Qualitatively, for each Q with < Q < 771,^ the illustrated dependence is the same. 
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Fig. 5 
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